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ABSTRACT 

We investigate the physical situation above the pulsar polar cap (PC) where 
the accelerating primaries (electrons) are not capable of producing sufficient num- 
bers of electron-positron pairs at low altitudes (within ~l-2 stellar radii above 
the PC surface) to screen the accelerating electric field, E\\, and continue acceler- 
ating up to, at least, the very high altitudes nearly approaching the light cylinder. 
We derive an analytic solution for the E\\ valid at high altitudes. The solution 
is based on the physical condition of asymptotic vanishing of the rotationally 
induced transverse electric field within the magnetic flux tube. This condition 
constrains the asymptotic value of the effective space charge that determines the 
distribution of electrostatic E\\ within the magnetic tube. Our estimates of low- 
to high-altitude values of E\\ imply the occurrence of a regime of primary accel- 
eration (with the characteristic Lorentz factor up to ~ 1 — 2 x 10^) all the way 
from the PC to the light cylinder limited by curvature-radiation reaction. In this 
model the primary outflow becomes asymptotically force-free, and may turn into 
a relativistic wind beyond the light cylinder. Such a solution will apply to both 
older pulsars producing only inverse Compton scattering pairs and younger very 
high B pulsars (magnetars). We suggest that pulsars, which are lying below the 
pair death line, may be radio-quiet 7— ray sources. 

Subject headings: acceleration of particles — gamma rays: theory — pulsars: 
general — radiation mechanisms: nonthermal — stars: neutron 

1. INTRODUCTION 

It is believed that one of the fundamental properties of isolated pulsars is the occurrence 
of electron-positron pairs (see Sturrock 1971, Ruderman & Sutherland 1975, and Arons & 
Scharlemann 1979) above their polar caps (PCs). Recent detailed investigation of the ac- 
celerating conditions in isolated neutron stars (NSs) implies, however, that the acceleration 



^Present address: ManTech International Corporation, Lexington Park, MD 20653 

^Laboratory of High Energy Astrophysics, NASA/Goddard Space Fhght Center, Greenbelt, MD 20771 



- 2 - 



of primary particles above pulsar PCs may not necessarily result in electron-positron pair 
formation (see e.g. Harding & Muslimov 1998, [HM98]; 2001, [HMOl] and 2002, [HM02] Hib- 
schmann & Arons 2001, and Harding, Muslimov & Zhang 2002 [HMZ02]). HM02 studied 
magnetic pair production by both curvature radiation (CR) and inverse Compton (ICS) pho- 
tons and the possible screening of the accelerating electric field by pairs. They showed that 
younger pulsars (with ages r < 10^ yr) are capable of producing sufficient numbers of pairs 
through CR to completely screen the electric field at low altitudes, while many older pulsars 
are able to produce only pairs through ICS in numbers that are insufficient for screening. 
For these older radio pulsars below the CR death line, as well as an unknown number of 
even older (and radio quiet) pulsars below the ICS death line, the primaries may keep accel- 
erating and emitting high-energy photons up to very high altitudes without any significant 
pair production. In the present study we develop the electrodynamic model which may be 
applicable to the description of high-altitude primary acceleration and 7-ray emission in 
both pair-starved radio pulsars and pairlcss, radio-quiet pulsars. We derive a simple ana- 
lytic solution corresponding to the steady-state flow of primaries and discuss the physical 
constraint on the solution near and beyond the light cylinder. This constraint determines 
the high-altitude solution that can be analytically matched with the known solution at low 
altitudes. This acceleration model can be used for older (r ~ 10^ yr) radio pulsars, including 
the majority of millisecond pulsars, as well as pulsars that are lying below the theoretical 
pair death lines (say in the P-P diagram) and may manifest themselves only as 7-ray sources 
(unidentified EGRET 7-ray sources?). Thus, the mode of operation of an isolated NS with 
dominated by high-energy emission and no pairs can be tested observationally. 

Some of the unidentified, and so far, radio quiet, 7-ray sources detected by the EGRET 
telescope of CGRO have been associated with the Gould Belt (Grenier 2000, Gehrels et al. 
2000). These faint sources have softer spectra and seem to be different from unidentified 
EGRET sources along the Galactic plane and from the known 7-ray pulsars. Harding & 
Zhang (2001,[HZ01]) explored the possibility that radio-loud 7-ray pulsars viewed at a large 
angle to the NS magnetic pole could contribute to this new population of Galactic uniden- 
tified EGRET sources. They showed that the bulk of the Gould-Belt 7-ray sources could 
be pulsars if their emission is seen at large angles to their magnetic axes, such that we are 
missing the radio emission and hard 7-ray emission at lower altitudes but detecting only the 
higher- altitude off-beam emission. Due to the flaring of the dipole fleld lines, this emission 
may be seen over a large solid angle, far exceeding that of the main beams. HZOl and Hard- 
ing et al. (2004) estimated that many of the EGRET sources associated with the Gould 
Belt may be such off-beam 7-ray pulsars. The population of intrinsically radio-quiet pulsars 
below the ICS pair death fine could also be contributing to the unidentified EGRET sources 
in the Gould Belt. Their 7-ray emission, while less luminous, would also be seen over a wide 
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range of angles and the acceleration model developed in this paper can be used to elucidate 
their predicted characteristics. We should point out (HZOl) that off-beam emission is not 
expected in outer gap models (see e.g. Yadigaroglu & Romani 1995, Cheng & Zhang 1998), 
which predict that most low Galactic latitude unidentified EGRET sources in the plane are 
radio-quiet, Geminga-like pulsars. 

An acceleration model for the younger pulsars above the CR pair death line has recently 
been investigated by Muslimov & Harding (2003, 2004, [MH04]). Even these pulsars that 
have large pair multiplicity and complete screening of the electric field still have a narrow 
slot gap along the last open field lines where no pairs arc produced and where the field 
is screened. High altitude particle acceleration is also possible for such pulsars, but MH04 
showed that the solution for the slot gap field differs significantly from the solution presented 
in this paper for high-altitude acceleration over the entire open field region above the PC. 

The paper is organized as follows. In §2 we introduce the electrodynamics of relativistic 
primaries accelerating above the polar cap of the NS. The constraint on high-altitude effective 
space charge is discussed in §3. In §4 we present the solution for the accelerating electric 
field valid at high altitudes and discuss the construction of general solution extending from 
the polar cap surface up to nearly the light cylinder. Finally, in §5 we discuss implication of 
our result for 7-ray pulsars, and draw our principal conclusions. 



2. The Main Electrodynamic Equations 

Let us consider the magnetosphere of a rotating NS and assume that in the frame of 
reference rigidly corotating with the NS the magnetic field is stationary. Then, the general 
relativistic Maxwell's equations yield (see MT92, MH97) 

E+ — (u- w) X B = --V$, (1) 
ac a 

where E and B are the electric and magnetic fields defined in Zero- Angular-Momentum- 
Observer (ZAMO) frame of reference (see Macdonald & Thorne 1982 for more detail); u = 
f2 X r and w are the velocity of the frame of reference rigidly corotating with a NS (Q 
is the stellar angular velocity), and the differential velocity of rotation of inertial frame 

of reference, respectively, relative to the distant observer; $ is the scalar potential; a = 
a/1 — s/f] {s = rg/ R, Tg is the gravitational radius of NS of radius R, and rj = r/R) is the 
so-called general-relativistic 'lapse function'; and operator V should be taken in curvilinear 
coordinates. Taking the divergence of eq. (1) and making use of Maxwell equation 



V • E = 47rp, 



(2) 
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we get the Poisson's equation for the scalar potential $ 



V. (iv*) =-47r(p-p^,), (3) 



a 



where 



— (u - w) X B 

a 



(4) 



is the general relativistic counterpart of the Goldreich- Julian (GJ) charge density (cf. Gol- 
dreich & Julian, 1969), and p is the actual charge density of electrons determined by their 
relativistic flow along the magnetic fleld lines and which is fixed by the condition E\\ = at 
the stellar surface (see MT92 for details). 

In this paper, in most of our calculations (unless otherwise stated) we assume that the 
NS's magnetosphere has two distinctive regions: the corotating region with closed field lines 
filled with charged plasma, and the open-field line region extending beyond the maximum 
radius of corotation (in cylindrical coordinates with the 2;-axis along the NS's rotation axis), 
the light-cylinder radius R\c ~ c/Q. As in our previous studies, we will be working in 
magnetic spherical polar coordinates {rj — r/ R,9, (f)pc), and denote by x the pulsar obliquity 
(angle between the NS rotation axis and magnetic dipole moment). We will also refer to 
a 'normal polarity' pulsar as one having 0° < x < 90° (north magnetic pole near north 
astrographic pole: f2 • m > 0, where m is the NS magnetic dipole moment), and a 'reversed 
polarity' pulsar as one having 90° < X ^ 180° (north magnetic pole near south astrographic 
pole: 17 • m < 0). 

In coordinates (77, 9, (ppc) eq.(3) takes the following form 

1 



where 



-Mp-Pgj): (5) 



is the angular part of the Laplacian. 

Note that the LHS of eq. (1) can be treated as the electric field seen by the observer 
who is rigidly corotating with the NS (who moves with the velocity u — w, i.e. with the 
magnetic field lines within the light cylinder), 

E' = E+ — (u- w) X B. (7) 

ac 

In our model calculation we shall assume that the magnetic field has pure dipole geome- 
try, at least well within the light cylinder. It is important that the static-dipole magnetic 
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field and the small-angle approximation (see below) become inaccurate for a rotating NS 
magnetosphere at very high altitudes as we approach the light cylinder. In addition, as 
we approach the light cylinder, wc will see (in ZAMO frame of reference) a strong trans- 
verse electric field induced by the rotation of the magnetic flux tube. Near and beyond the 
light cylinder the magnetic flux tube should rather rotate at a retarded u < c, sweeping 
back in the direction opposite to rotation, thus partially reducing the electric field induced 
by rotation. Also, this electric field may be reduced by the effect of particle drift across 
the magnetic fiux tube, but this effect is less pronounced for the entire open field line re- 
gion compared to e.g. the narrow slot-gap region (MH04), where this drift may dominate. 
Note that, at the very fight cylinder, where the velocity of the magnetic tube is formally 
approaching, and becomes equal to, the velocity of light, the effective transverse electric 
field as seen by the observer corotating within the magnetic tube E_|_ + c"^(u — w) x B 
should vanish to ensure that E' is finite, because E'_|_ r[E_|_ + c~^{u — w) x B], where 
r ^ [1 — (■u/c)^]""^''^ oo as u c. More generally, we suggest that the condition of vanish- 
ing of E_L + c~^(u — w) X B is the basic physical condition (see also MH04) for sustaining the 
continuous steady-state fiow within the magnetic fiux tube and which also ensures that the 
last open field fine surface will remain equipotential (MH04). Note that p (as well as all other 
electro dynamic quantities) is defined in ZAMO frame of reference, and at higher altitudes it 
is mostly determined by the rotation of the magnetic flux tube. Asymptotically, at higher 
altitudes, the ZAMO sees the electric field dominated by the (transverse) rotation-induced 
component E ^ — (u/c) x B, which is associated with nearly GJ (fictitious) charge density, 
p = (l/47r)V ■ E f» — (l/47rc)V ■ (u x B) = p^j. Thus, in a steady state, within the magnetic 
flux tube, the transverse electric field induced by rotation of the tube vanishes. 

Ex + - X B 0, (8) 
c 

and the electric field within the magnetic tube is essentially determined by the distribution 
of the actual space charge within the tube. Here, for brevity, we have assumed that u':$> w at 
high altitudes. The condition (8) allows us to derive the high-altitude solution for without 
explicit calculation of p (and also B) at very high altitudes. In contrast, at very low altitudes 
the rotation only slightly perturbs the steady-state solution and may result in establishing 
a weak current along the boundary field lines without violating the equipotentiality of the 
boundary (sec discussion followed by eq. [35] in MH04). In this case, p is simply determined 
by at the stellar surface, and the electric field, E, is mostly determined by the distribution 
of physical charges ejected from the PC surface. Note that condition (8) is called the 'ideal 
MHD condition' (see e.g. Mestel, 1999), though here it simply refiects the induction of an 
electric field by a rotating magnet and may be valid even for a single test particle moving 
along the field fine. [Strictly speaking, the MHD approximation is only valid in coUisional 
plasmas and also may fairly well describe the perpendicular (but not the parallel!) motions 



6 



of coUisionless plasmas.] 

Above the PC surface (and up to very high altitudes) the radiation drag (e.g. due to 
curvature radiation) may become important, and the motion of primaries along the magnetic 
field lines can ignore inertia and be satisfactorily approximated by the balance of eE\\ and 
the drag force. In the end of Section 4 we use this approximation (see also MH04) to 
estimate the characteristic Lorentz factor of a primary moving in the regime limited by the 
curvature-radiation reaction. 



3. The Constraint on Effective Space Charge at High Altitudes 

The condition (8) implies that at very high altitudes, near the light cylinder, maintaining 
a steady-state and continuous flow of charged particles along the magnetic stream lines 
requires that the effective space charge density, |Ap| = \p^j — p|, remains -C We 
assume therefore that the electrostatic solution which will be valid up to high altitudes 
within the light cylinder should comply with the asymptotic condition (8). In this Section 
we shall illustrate that the high-altitude expression for the effective space charge density 
can be determined by employing the standard expression for p valid within relatively low 
altitudes (or, for the static magnetic field) and p^j calculated for pure dipole magnetic field 
rigidly co-rotating with the star. The only constraint we need is the asymptotic vanishing of 
E_L -I- u X B/c along the magnetic stream fines. It proves instructive to use the effective fiux 
of charges (both actual and GJ) through the cross-sectional area of a magnetic fiux tube (of 
dipole field) emanating from the PC, defined as (see MH04) 



where f{rf) is the general-relativistic correction factor (see Appendix A). 

The fiux of charges streaming with relativistic velocity through S{r]) can be written as 



where p is the local charge density (as seen by ZAMO). Then, using the explicit expressions 
for p (valid for low altitudes/slow rotation and coinciding with the actual charge density 
seen in the frame of reference corotating with the magnetic fiux tube) and p^j, for the fiux 
of actual and GJ charges we get the approximate (for illustration) expressions (see e.g. eqs 



5(77) = ^5(1) 



M,3 



(9) 



F = a{r])cS{rj)p{ri), 



(10) 



[20], [21] in MH04) 




(11) 
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and 




K 3 7] ^r,\ 3 nr. . 



(12) 



respectively. Here the superscript 'd' refers to the dipole field; Fq = —pQS{ri = l)c, po = 
VLBq/2-kc and n is the magnitude of the general- relativistic effect of frame dragging at the 
stellar surface measured in stellar angular velocity Q,. For most more or less realistic NS 
equations of state, k ^ 0.15 I45/RI; where I45 — I/IO^^ g-cm^, Rq — cm, I is the 

moment of inertia of NS of radius R) . 

In expressions (11) and (12) we omitted the general-relativistic red-shift corrections (factors 
H and / which are presented in Appendix A) and factor ^/l — 3^'^rj/4:7]\c in the expression 
for p that takes into account the weakening of the dipole magnetic field strength towards the 
magnetic equator (see eq. [21] in MH04). 

At the very high altitudes, nearly approaching the light cylinder, the magnetic stream 
lines begin deviating from those of a vacuum magnetic dipole, and, more importantly, the 
terms associated with the rotation become the leading terms in the expression for F (which 
will be different from formula [11] based on the static magnetic dipole approximation), so 
that F to comply with condition (8). We shall now demonstrate this on a simple 

model. We assume that near the light cylinder the open magnetic field lines can be roughly 
approximated as radial lines. This assumption seems to be reasonable and may correspond 
e.g. to the situation where beyond the light cylinder the kinetic energy density of particles 
(primary electrons) exceeds the magnetic energy density, so that the magnetic field lines 
become governed by the outflowing particles. Or, there may be a regime of force-free current 
discussed by Contopoulos, Kazanas, & Fendt (1999), and Mestel (1999) who neglected the 
particle inertia and assumed that the local current density is nearly equal to the GJ current 
density. As a result, in both cases, the magnetic field lines (in open field lines domain) of 
an initially dipole magnetic field are expected to straighten out near (or even well within) 
the light cylinder to become quasi-radial and then become dominated by the azimuthal 
component further out beyond the light cylinder. For the radial magnetic field, B = B^-e,. 
(in spherical coordinates r, 9, 4>pc, and with unit vector in radial direction er) we can get the 
following expression for the flux of GJ charges (here, for the sake of illustration, we ignore 
correction factors attributed to the static general-relativistic effects) 



In derivation of formula (13) we assumed that the open field lines of the poloidal magnetic 
field become nearly radial at 77 ~ 773, and that 



F: 



irad 



Fq (cos X COS 9 + sin x sin 9 cos 0pc) . 



(13) 



GJ 




(14) 
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where Bq is the surface strength of dipole magnetic field at the magnetic pole. Note that 
eq. (13) implies that the magnetic stream lines are radial lines {6 = const), whereas eq. 
(12) implies that, on magnetic stream lines (see formula [19] below), 9 = sm~^{^^yr|/r|lc). 
From expressions (12) and (13) one can see that F^^ can satisfactorily match the F^^^ for 
rj < (l/3)?7ic and 9 < 1/3. Due to the charge continuity equation, the quantity p/B (for static 
magnetic field) should be constant along the magnetic stream lines. This means that for 
the static magnetic field with dominating radial component, p oc S oc r]^"^, and therefore 
F'^^ (oc prj'^) should also be constant along the magnetic stream line. Even though this 
result is trivial, it illustrates the general fact that, in a steady state, the fiux of charges 
streaming within the magnetic tube remains constant and does not depend on the geometry 
of the tube. To derive a general form of F at all altitudes which complies with condition 
(8), we match F'^ with F^^'^ at high altitudes. In fact, this very general reasoning allows 
us to restore the high-altitude flux F seen by ZAMO. For example, for F we can write the 
following expression containing the same radial dependence as in F^^ at rj ^ rjc (where rjc is 
some adjustment radial coordinate to be determined), 

3 

Vc) - 7^[{V - Vc)h{r] - Vc) + M^'^ COSX+ 



Fo 

3 
2 




2r7ic 



sm X cos 0pc 



(15) 



where, for simplicity, we used step function h{x) {h{x) = 1 if x > 0; 0, otherwise), so that 
for 77 < we arrive at formula (11), while for 77 > 77c we get the asymptotic expression for F. 
Note that for our purpose the specific functional form of adjustment of to its asymptotic 
expression at very high altitudes is not important, because the only physical quantity we 
actually need is the effective flux, F^_^ — F, within the magnetic tube. Since the effective 
flux is constant at higher altitudes, it can be determined just at very high altitudes where 
asymptotic condition (8) is satisfied. Thus, using expressions (12) and (15) we can get 



AF 



F = Fn 



- -T - 3 ^^(^7 - ^c) 



2r7ic 
3 



77' 

cos x+ 



For T) < Tie formula (16) reduces to 
AF = Fo 
3 



[77 - 1 - (77 - 77c)/l(77 - 77c)]C^ 

[07 - 1 - (07 - 07;)/i(77 - 77c)]Csinxcos0pc| . 



(16) 



1 - ri)e 



cosx 



^o^iVv- l)sinxcos( 



(17) 
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and for rj > rjc we get 



AF 



3 



K 1 - 



COS x+ 



Oo{^/rk- l)Csinxcos( 



'pc ( , 



(18) 



where the general- relativistic factors /(I) and H{1) are temporarily omitted for the sake of 
simplicity. 

Formulae (17) and (18) turn out to be very general and illustrate that AF gradually varies 
from low- (eq. [17]) to high- altitudes (eq. [18]), and that at very high altitudes AF is fully 
determined by F^^ — F'^ fixed at low altitudes and perfectly conforms to condition (8). Also, 
formula (18) implies that AF should saturate with altitude (see also MH04), and that this 
saturation occurs well within the light cylinder. Using the high-altitude expression for AF 
given by (18) we can now find the appropriate solution for $ which will depend on the 
adjustment parameter, radial coordinate rjc- Then, by matching the low-altitude solution for 
$ with the corresponding high-altitude solution, we may determine the parameter rjc (see also 
MH04). We must note that for estimating e.g. the energetics of acceleration of primaries 
we need to know only the low- and high- altitude solutions, since [$| is a monotonically 
increasing function of r] along the magnetic stream lines. 



4. The High-Altitude Solution for Potential $ Above the PC 

We begin with the general formula for the static dipolar magnetic field line (which is 
formally valid up to arbitrary high altitudes within the light cylinder), 

where = [^-R/c/(l)]^''^, ?7ic = c/QR, and ^ is the dimensionless magnetic colatitude of 
open field lines (^ = 1 corresponds to the last open field lines, and = corresponds to the 
magnetic axis). [In a small-angle approximation, which implies that r]^^/r)\c -C 1, we can 
write 9 a; ^^/rjJr|^c■] 

At high altitudes, but still within the light cylinder, where 3> d"^ /drf, we can use 
expression (18) and relationship p^j — p — AF/acS{rj) to rewrite eq. (5) as follows 

Vn$ = — (^cosx + 'Bsinxcos^pe), (20) 
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where J^q = — 2$o//(l), = {Q.R/c)Bq R. We shall introduce new dependent variables, 
yi = ?7$i and 1/2 = Tj^2) and independent variable x = cos6 = ^1 — a'^^'^ (where a = \/rjJr]\c 
and approximation sin(^^o) ~ C^o is used). Then, the expressions for A and B can be written 
in general form as 

A = Ai + (21) 

B = BiVT^ + B2xVl^. (22) 

We will search for the solution of eq. (20) along the magnetic stream lines as having the 
form 

$(?7,^,0pc) = JP'o($icosx + $2sinxcos0pc), (23) 
Now eq. (20) translates into 



d_ 
dx 



:i-x' 



d_ 
dx 

^dy2 
dx 



'l-x' 



dx 



Ai + A2x'^, 



y-2 



1-^2 



BiVl^ + B2xVl^ 



X^ 



(24) 



(25) 



The general solutions of eqs (24), (25) vanishing at the last open field line (equipotentiality 
of the last open field lines, = 1] = 0) read 



yi = - ( ^1 + - A ) In 



l+x 

1 + Xq 



\^A2{i-e), 



y2 = oa/— C 

2 V ^Ic 



Bi 



1 +a:o 
l + x 



l]+-B2{xo 



l + x 



X 



(26) 
(27) 



where Xq — Vl — = ^/^ic- 



In a small-angle approximation [j]^ jrixc -C 1), when a; ~ 1 — [77/27710]^^, expressions (26) 



and (27) reduce to 



Z/2 



1 f rj 



3/2 



{B,+B2M^-e)- 



8 V^ic 

The general solution for $ valid at high altitudes within the light cylinder, reads 



(28) 
(29) 



2 Ml + -^2 



Vic 



V 



In 



l + x 

1 + Xo 



+ 



^A2{i-e) 



Vic 



COSX+ — 



1/2 



-Bi 1- 



1 + 3^0 

l + x 



+ -B2 { X - Xo 



1 + a;o 
l + x 



sm X cos 0pc 



(30) 
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In a small-angle approximation the solution (30) reduces to 

1 



-$o^^[(A + ^2)cosx+ 



^^0 VvfWW)il3i + 62)^ sin X cos . 



'pc 



(31) 



Using expression (18) we can now identify the explicit values of Ai, A2, Bi, and B2 con- 
forming to the condition (8) and that determine the high-altitude solution for $ (see also 
MH04), 

(32) 



VcJ ^ /(I) V 

^i^?(v^-i)4^^, 



3,_, ,ii'(l)l 



!^ 0. 



(33) 



Here we restored general-relativistic factors /(I) and -f^(l) that were omitted in all corre- 
sponding expressions in previous Section used for illustration. In Appendix B we present 
general and approximate expressions for the functions determining the radial dependence of 
$ and £'11 . 

In a small-angle approximation (this is good for small %, ^ and at high altitude) we 
arrive at (cf. eq. [31]) 



1 fnny ^ 

1/2 



K\l 5- I COS X-F 



3 {VLR\ ' 
-(— j A(^- l)^sinxcos0pc 



(1-a- 



(34) 



The corresponding expression for E\\ — —{l/R)d^/dr] reads 



E\\ — -E'||,high 

1 



3 fnuy Bo 



16 V c ; /(I) 



A(l + 2e')esinxcos0. 



pc 



(1-a- 



(35) 



In expressions (34) and (35) X ^ H{1)/ ^/f{^ ^ 1- 0.635-0. 16£2-0.5fi;(l-0.63£-0.2l£2). 
By comparing the low-altitude solution (see eq. [37] below and e.g. eq. [13] in HM98) and 
approximate high-altitude solution (34) one can see that the potential $ practically saturates 
at 77 > 77c, and the small growth of $ occurs over a length scale of order of the light cylinder 
radius and is determined by terms oc rj/rjic- 
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Note that the radial dependence in formula (30) is mostly determined by the effect of flaring 
of the magnetic field lines (for large ^), so that near the very light cylinder (on the last 
open field lines), where this effect is especially profound our analytic solution becomes less 
accurate. However, for most pulsar parameters our solution is quite satisfactory within 
~ 0.5 — 0.7 of the light cylinder radius, where approximate formulae (34) and (35) can be 
used. 

We can now calculate electron acceleration up to very high altitudes by employing the 
following simple formula for E\\ 



where i?||,high is determined by expression (35), and £^||,iow is given by (see e.g. formula [13] 
in HMQsj 



where /3 ^ l + 0.63x + 0A6x'^ -l.Ux^ + {3K/ri^-^){l -0.83x -OMx^ + 0.61x^), and a; = e/rj. 

In formula (36) the parameter rjc can be determined as the best-fit parameter that pro- 
vides a smooth transition between the low- and high-altitude values of E\\. As we have 
discussed in our previous study (sec MH04), construction of a general solution for E\\ differs 
markedly for the regions with the "favorably" and "unfavorably" curved field lines introduced 
by Arons & Scharlemann (1979). The region with favorably curved field lines corresponds 
to that with cos^pc > (electron acceleration for normal polarity pulsars, and positron 
acceleration for reversed polarity pulsars) , and the region with unfavorably curved field lines 
corresponds to that with cos0pc < (positron acceleration for normal and electron accel- 
eration for reversed polarity pulsars, respectively). In the classical electro dynamic model 
of Arons & Scharlemann (1979) the accelerating electric field above the pulsar PC is de- 
termined by the component oc sinxcos0pc, and therefore the electrons or positive charges 
(positrons) may accelerate in their model only along favorably or unfavorably curved field 
lines, respectively. In our model the presence of the component oc cos x (determined by the 
general-relativistic frame-dragging effect) tends to make all field lines favorably curved in the 
sense that the charges of the same sign are accelerated outward over the PC. However, for 
very large obliquities and/or very short spin periods the field lines may become unfavorable 
for electron acceleration, because of the negative contribution from the term oc sin x cos 0pc 
(cos0pc < 0). Earher (MH04) we have mentioned that for nearly ahgned (0°<x<60°) or 



£;||,iowexp[-(r7 - l)/{r}c - 1)] + -Ey^high, 



(36) 




(37) 
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nearly anti-aligned (120° <x < 180°) cases, there will be continuous acceleration of electrons 
from the PC surface up to very high altitudes on both favorably (cos 0pc > 0) and unfavor- 
ably (cos0pc < 0) curved field lines. For nearly orthogonal rotators (80°<x~100°) there 
will be continuous acceleration of electrons/positrons from the PC surface up to very high 
altitudes on favorably curved field lines for normal/reversed polarity, and continuous accel- 
eration of positrons/electrons on unfavorably curved field lines for normal/reversed polarity. 
At intermediate inclination (60° < x < 80° and 100° < x < 120°), there will be continuous accel- 
eration of electrons on favor ably /unfavorably curved field lines for normal/reversed polarity. 
Let us now summarize the characteristic values of rjc in formula (36) for the situations where 
in this formula either the term oc cosx or oc sinxcos0pc dominates. We will present our 
estimates for the value of ^ ~ 0.5. For the case where the term oc cosx formula (36) 
dominates, r/c ~ 3 — 4 and rjc ~ 10 — 15 for the millisecond pulsars and relatively long-period 
(P ~ 0.5 — Is) pulsars, respectively. For the case where term oc sinxcos^pc dominates, 
r^c ~ 3 — 4 (millisecond pulsars) and rjc ^ 15 — 20 (pulsars with P ~ 0.5 — 1 s). Thus, for a 
given spin period the parameter rjc is, generally, a function of x and (ppc- 

Let us now estimate the characteristic electron Lorentz factor, 7, in the regime of acceleration 
above the pulsar PC limited by the CR reaction, i.e. where the condition 



2 



e|^|||-^^7^ (38) 

is satisfied. Here Pc ~ (4/3) (-R/C) ijlVic)^^'^ is the characteristic curvature radius of the field 
lines (0 < ^ < 1). 

By inserting expression (35) for E\\ into expression (38), we obtain 



A 



7^2xlO^(^J p^'jr ] Kx.iA,.y'\ (39) 



(^l--^j (l + f)cosx+2(v^-l) y—j -(l + 2aesinxcos0pe 



(40) 

where B12 = Bq/W^ G, «;o.i5 = «;/0.15, Pq.i = i'/O.l s. For the ms (e.g. P = 3 ms, 
B12 — 3 • 10~^, and Rq — 1) and for 1 s (P12 — 1,Rq — 1) pulsars we get 7 ~ 10^. 

The characteristic energy of CR spectrum, £cr, is given by 

S^^/m^c^ = 3Ac7V2 Pc, (41) 
where Ac = 3.9 x 10~^^ cm is the electron Compton wavelength. By inserting expression (39) 
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into formula (41), we get 

fTc. ^ 400 MeV ( ^) { ^'^^^^A ' aV\ (42) 

Thus, for the ms and ~ 1 s pulsars Ecr ~ 600 and ~ 10 MeV, respectively. Thus, the off-beam 
spectrum of CR for the ms pulsars occurs at EGRET energies, but may be at a flux too 
small to detect except in the case of very nearby sources. The on-beam CR emission from 

ms pulsars which is radiated nearer the PC (Harding, Usov & Muslimov 2004) occurs at 1-10 
GeV. since i^mow ^ ^||,high- It is interesting to point out that near the light cylinder, where 
the radius of curvature of open field lines may substantially increase, the curvature radiation 
may cease. However the shape of the curvature-radiation spectrum and peak spectral energy 
will be determined by the harder photons produced at low altitudes. 

Let us also estimate (within the small-angle approximation) the characteristic luminosity 
of high-energy off-beam emission. For example, for cosx ~ 1, we can write 

^HE = /he ("^c^) c (^er^) / rf0pc / sin Ode 

Jo Jo 

« (mc^) c {n^T]') {2nrl,) [\mdi, (43) 

^0 

where /y^i (< 1) is the efficiency of conversion of electron kinetic energy into high-energy 
(curvature) emission (which may be calculated by integrating the CR power), 9aiax — C^oy/v^ 
T^pc — (^oR is the PC radius, and Ue ~ QSo(l — K)/{2T:cerf) is the electron number density. 

By substituting expressions (39) and (40) into formula (43), and using kq.is = 1, we get 

rj5/4 pl5/4 / \ 1/4 

L,,«2xlO-erg.s-/.,i^^0FUr " ^^^^ 

Note that, according to formula (44), L-^^ cx L^J^ , where L^d (~ 5 ■ 10'^'^ erg - s"^ -^6 -^oTi^) 
is the pulsar spin-down luminosity. We can also estimate the efficiency of conversion of 
pulsar spin-down luminosity into high-energy luminosity, 

/„, = ^ = 0.07 (A £.tS. (45) 

sd \ / 

The estimate (44) is roughly consistent with the expected luminosity of off-beam pulsars (see 
HZOl). However, for the quantitative modeling of high-energy fluxes from these hypothetical 
sources the above formulae, based on a small- angle approximation, may not be accurate and 
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the appropriate detailed treatment will be needed. According to (45), the efficiency of high- 
energy emission significantly increases for long-period pulsars (both with regular, ~ 10^^ G, 
surface magnetic field and magnetars) and also for the very low-field ( < 10* G) ms pulsars. 
Note that in estimates (39), (42), (44) and (45), according to a small-angle approximation, 
1] <^ ^ic/C^ 7^ 0). This means that for a nearly aligned pulsar and sufficiently small 
magnetic colatitudes (say ^ ~ 0.1 — 0.2) these estimates are formally valid for very large 
radial distances, rj^rjic- It is interesting that for small ^ the values of 7 and (see [39] and 
[42]) are proportional to 1/^ which implies the substantial hardening of the CR spectrum 
towards the magnetic axis. Obviously, ^ cannot be arbitrary small, with the lower limit set 
by condition Ljj^ < Lg^- 

In this paper we discuss the possibihty of occurrence of the regime of primary accelera- 
tion up to high altitudes without any significant pair formation at low altitudes. However, 
the pairless regime of acceleration of primaries may not be stable (at least, for some mag- 
netic colatitudes) against ^ + X —>■ e'^ + e~, 7 + 7 + e~, or even (for magnetars) 
^ + B — > e'^ + e~ processes. We plan to address this issue and calculate the emergent 
high-energy spectrum for the relevant scenarios separately. We expect, therefore, that the 
high-altitude pair formation, if it occurs, may result in effective narrowing of the accelera- 
tion region and associated weakening of the accelerating electric field. It is unlikely that the 
regime of pairless acceleration can be completely shut down by high-altitude pairs, though. 



5. Discussion and Conclusions 

We have considered the possibility of continuous acceleration of primaries (electrons) 
on open field lines above the pulsar PC up to very high altitudes, nearly approaching the 
light cylinder. Such a possibility may occur in pulsars which lie below the CR pair death 
fines (see e.g. HM02, and HMZ02), because the efficiency of pair formation by ICS may 
not be sufficient to screen the accelerating electric field. These results will apply to old (age 
~ 10^ yr) and ms radio pulsars, as well as to radio-quiet pulsars that are not producing 
any pairs. Well within the light cylinder the magnetic field is well described by pure dipole 
structure, and primary fiow is governed by the magnetic field. However, as we approach 
the light cylinder, the particle stream lines may drastically deviate from the magnetic field 
lines of a static dipole, and beyond the light cylinder the fiow of relativistic particles may 
govern the magnetic field. As a result, the distribution of GJ and actual charge densities 
near and beyond the fight cylinder may significantly differ from that for low altitudes. By 
considering the magnetic flux tube of a rotating NS and taking into account the fact that at 
high altitudes, the strong transverse electric held as seen by a ZAMO and which is induced 
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by the rotation of the magnetic flux tube, vanishes within the tube, at very high altitudes 

the effective charge density |Apeff| = IPgj ~ p\ IPgjI- ^^^^ condition which provides 
electrostatics within the magnetic flux tube, and allows us to extend the low-altitude solution 
to much higher altitudes nearly approaching the light cylinder. Thus, we derived explicit 
expressions for the scalar potential and accelerating electric fleld valid through the very high 
altitudes, where the often- used small-angle approximation (condition i^rj/rj^^ -C 1) becomes 
invalid (for not very small ^), and the corresponding low-altitude solution is not applicable. 

We have shown that the low-altitude solution for the scalar potential $(?]) makes a 
smooth transition to the high-altitude solution constructed in the present study for which 
the corresponding charge density of relativistic primaries may match the wind solution akin 
to that derived in the force-free case (Contopoulos, Kazanas, & Fendt 1999), thus connect- 
ing the near-surface pulsar electrodynamics to a global solution for the magnetosphere and 
wind. The high-altitude solution presented in this paper allows us to estimate energetics 
of primaries that are not capable of producing pairs at very low altitudes above the PC 
and therefore accelerate up to the very high altitudes. For example, wc have demonstrated 
that these primaries may achieve signiflcant Lorentz factors (up to a few times 10^) and 
emit 7-rays in the regime of acceleration limited by the CR reaction. The characteristic 
length scale of high-altitude acceleration is roughly equal to the light cylinder radius, and 
the high-energy luminosity scales as L^^^. 

In conclusion, we should point out that most of the difficulty in constructing the quan- 
titative model for high-altitude primary acceleration arises from the lack of the reliable 
knowledge about the magnetic fleld structure of rotating NS near and beyond the light 
cyhnder. We have accentuated in this study that, at high altitudes, the rotation of mag- 
netic flux tube fully determines the distribution of the electric charge density within the 
tube as seen by the ZAMO (or, in a flat space-time, by the observer in the 'lab' frame 
of reference), and demonstrated the electrostatics within the magnetic flux tube which is 
important for the calculating the high-altitude accelerating electric fleld. Such a regime of 
steady-state acceleration of primaries is crucial for constructing a self-consistent steady-state 
global magnetospheric model. 

We acknowledge support from the NASA Astrophysics Theory Program through the 
Universities Space Research Association. 
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A. General-relativistic Corrections 

In this Appendix wc present explicit expressions for the physical components of the 
dipolar magnetic field around NS (in Schwarzschild metric). Also, we present some important 
general-relativistic factors entering the main electrodynamic equations in curved space-time. 
We employ the magnetic spherical coordinates r, 9, and 0pc, and use the standard notations 
for the fundamental physical constants. Other notations and shorthands are: 

R— the NS radius 

Tg— the stellar gravitational radius 

77 = r/R— the dimensionless radial coordinate 

e — Tgl R— the NS compactness parameter 

X — e/r] 

a = v^l — X— the general-relativistic 'lapse function' 

K ^ 0. 1 5 /45/-Rg — the frame-dragging effect at the NS surface measured in fl, angular velocity 

of stellar rotation 

J45 = J/IO^^ g ■ cm^ 

/— moment of inertia of the NS 

Re = R/W^ cm 

The physical r— , 9— components of dipolar magnetic field are 

Bf^Bo^^cos9 (Al) 



and 



-2M + 



^ sin^, (A2) 



jy3 



respectively, where Bq is magnetic field strength at the magnetic pole. Here 

fiv) = -^M^-x) +x{l + x/2)]. (A3) 

For large r] (small x) one can use approximate formula 

f{rj) ^ 1 + 0.75x + 0.6x^ (A4) 

The correction factor H{r]) entering the general-relativistic expression for the Goldreich- 
Julian charge density reads 

Hi,) = X - .[x/ef + 1 [,_l^^f{,;^ (A5) 

The approximate expression (for small x), 

H{ri) ^ 1 - 0.25a; - 0.16x^ - 0.5K{x/ef{l - 0.25a; - 0.21x'^). (A6) 
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B. Radial Dependence of $ and E\\ 

By incorporating expressions (32)- (33) into eq. (30) we get 

$ = i$o^o(-^i cos x + sin x cos 0pc), (Bl) 

where 

^,=4(^, + lA)^l.>(i±^)+HA(l-«. (B2) 

and 

^= - (^) ^' 

where x = -^1 — ri^'^/rjic, and Xq = x(^ = 1). 

Here we present expUcit formulae for the functions JF^ and JF2 defined above, and their 
derivatives with respect to rj. Also, we express J-'i and J^2 as a power series in ry/r/ic, up to 
the quadratic terms. We use the same notations as in Appendix A, along with the variables 
Vic — c/D,R— the dimensionless (cylindrical) radial coordinate of the light cylinder 
^ = 9{r] — 1)1 9q— the dimensionless magnetic colatitude (see eq. [19]) 



The 77-derivatives of J^i and determine 



--^01 

2 R ° 



9^1 ^0^2 . , 

— — COSX+ ^— smxcos0pc . , 
orj or) ' 



and can be expressed as 
9 In 77 



2 ( ^1 + -^2 
A (4^ In 



xo{l + xo) x{l + x) 77 l + Xo 



+ 



1 + x 



V 1 + xo 




2{xo -x) - 



V_ f 1 + 3^0 ^2 _ }_ 
\x{l+x) Xq 



(B4) 



(B5) 



91n?7 y 1] 1 + X 

For Tj/rjic'^ 1/3 the expressions (B2)-(B3) can be approximated by the following formulae 

2 



(B6) 



^1 



8771c 



24 V?7i, 



+ 



3nRH{l) 

8V7(Ty' 



1-^c) 



i+e' + ^-(i+f +a 
9771c' 



(B7) 
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1 + 1^(1.20 



8 V^ic 



1 + + 



(B8) 



977 ~ 8 c 



5Qi?//(l) 



(B9) 



dri 



3 

16 



77(1) 



X 



?7ic 



(BIO) 
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